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The dc-magnetization of the unique S=l/2 kagome antiferromagnet Herbertsmithite has been 
measured down to 0.1 K. No sign of spin freezing is observed in agreement with former fiSR and ac- 
susceptibility results. The low temperature magnetic response is dominated by a defect contribution 
which exhibits a new energy scale ~ 1 K, likely reflecting the coupling of the defects. The defect 
component is saturated at low temperature by > 8 T applied magnetic fields which enables us to 
estimate an upper bound for the non saturated intrinsic kagome susceptibility at T = 1.7 K. 

PACS numbers: 75.30.Cr, 75.30.Hx, 75.50.Lk 



In triangular lattices, the frustration of antiferromag- 
nctic interactions associated to the enhancement of quan- 
tum fluctuations for S=l/2 spins was acknowledged long 
ago as a keypoint to stabilize novel ground states of mag- 
netic matter^. Numerous theoretical studies have since 
then emphasized the S=l/2 nearest-neighbor Heisenberg 
antiferromagnet on the kagome lattice (KAH), a net- 
work of corner sharing triangles. Although numerical 
approaches are complicated by the huge degeneracy of 
the system, it is believed that the ground state could 
be a unique realization of a disordered two dimensional 
quantum liquid at T = 0, with a surprisingly small gap, 
if any, to unconventional unconfined spinon excitations 
and a gapless continuum of non magnetic excitations^ii^ii^. 
Concurrently, growing efforts were made to identify a 
model frustrated compound and find evidences for such 
an exotic spin liquid ground state. Key features have 
emerged from these experimental investigations like the 
suppression of magnetic order at the energy scale of the 
antiferromagnetic interaction, the persistence of spin dy- 
namics at very low temperatures^ or a large density 
of low energy non magnetic states^. However, the frus- 
trated compounds studied so far show strong deviations 
from the ideal KAH (S> 1/2 spins, dilution of the mag- 
netic lattice, anisotropic interactions). Besides, they of- 
ten present marginal low T order or spin glass like be- 
havior which forbid a close comparison to theoretical ex- 
pectations*. 

Only very recently, Herbersmithite, ZnCu3(OH)6Cl2, 
a structurally perfect kagome antiferromagnet decorated 
by Cu^+ S=l/2 spins could be synthesized^. It belongs to 
a large compound family Zn^Cu4_2:(0H)gCl2 where the 
parent structure, clinoatacamite (x=0), is that of a dis- 
torted S=l/2 pyrochlore. The substitution of Zn^+ ions 
preferentially on the less Jahn- Teller distorted Cu^+ site 
located in between the kagome planes restores the three 
fold symmetry of the lattice for x > 1/3. Eventually the 
Herbertsmithite compound {x — 1) presents decoupled 
S=l/2 perfect kagome planes. Also, the magnetic order 
which sets in clinoatacamite at 19 K gradually disappears 



as X —^ 1 and for x = I, muon spin resonance (/iSR) in- 
vestigation^" has demonstrated the absence of any spin 
freezing at least down to 50 mK, an energy scale 4000 
times smaller than the main antiferromagnetic interac- 
tion (J ~ 190 K). 

Once Herbertsmithite is acknowledged to be the first 
good candidate for the realization of the KAH model, 
the magnetic susceptibility and heat capacity are the first 
quantities of interest as they straightforwardly probe the 
nature of the ground state, either magnetic or not, and 
the excitation spectrum. At low T, these thermodynamic 
quantities show respectively a Curie-like tai l^^d^ and a 
Schottky-typeiiii^ anomaly. It was soon recognized that 
Dzyaloshinsky-Moriya interactions, which are allowed in 
Herbertsmithite structure, can yield such a drastic in- 
crease of the low T susceptibilitjii^. However magnetic 
defects could also account for these features, a scenario 
sustained by recent NMR data— and neutron diffraction 
refinements of the structur o^^d^ . Despite a poor sensitiv- 
ity, these latter point at a large (6 — 10%) Cu/Zn inter- 
site mixing. This chemical disorder would likely reffect 
the finite energy of the Jahn- Teller process that selects 
the Zn substitution site. Both the resulting dilution of 
the kagome magnetic network and the interplane Cu^+ 
ions may contribute to the defect component. In this pa- 
per, we report on a detailed investigation of the Herbert- 
smithite magnetization at low temperature (T > 0.1 K) 
and up to moderately high fields {H < 14 T) . Both sets of 
data are consistently analyzed in terms of a large defect 
contribution. We show that the low T intrinsic suscepti- 
bility can be nonetheless estimated. 

A ZnCu3(OH)6Cl2 powder sample was prepared by the 
hydrothermal method described in Refi^iiS. Low tem- 
perature dc-magnetization (0.1 K< T < 3 K) was mea- 
sured in a home made SQUID magnetometer for fixed 
external magnetic fields up to 0.7 T. For T > 0.3 K, 
each data point was obtained by extracting the sample 
through the pick-up coils. To avoid heating effect, for 
0.1 K< T < 0.4 K, the sample was kept at a fixed po- 
sition in the pick-up coils and we measured the SQUID 
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FIG. 1: Molar dc-susceptibility of Herbertsmithite versus 
temperature in 0.1 T external field on a log-log plot. The 
local susceptibility measured by ^^Cl NMR line shift falls in 
the shaded areaiS. The dashed line is a Curie Weiss fit for 
1.5 K< T < 3 K (see text). Inset: T-dependence of the in- 
verse of the susceptibility at low temperature. 
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FIG. 2; (Color online) M/ H as a function of temperature for 
different H. Inset: M normalized by the saturated magneti- 
zation of one mole of S=l/2 spins (Cu'^"'') Msat = 5583 emu, 
at 0.2 K as a function of H . The data are extracted from the 
temperature scans of the main panel. The straight line is a 
plot of x(0.2 K)H. 



voltage variation as a result of the T-dependent sample 
magnetization. Besides, dc-magnetization curves were 
recordered versus field (0 < if < 14 T) at constant tem- 
perature in a commercial vibrating sample magnetome- 
ter (VSM). Standard SQUID data up to 5 T and for 
T >1.8 K were also used to complement and calibrate 
the low T data. 

The temperature dependence of the dc magnetic sus- 
ceptibility X of Herbertsmithite measured in a 0.1 T ap- 
plied field is presented on Fig. [1] in the whole studied 
temperature range (0.1 K< T < 300 K). At high tem- 
perature T >150 K), the susceptibility shows a Curie- 
Weiss behavior which yields the exchange constant J ~ 
190 Kiiiii. At lower temperature, the susceptibility in- 
creases much more rapidly down to ~ 0.5 K where it 
eventually flattens. Down to the lowest temperature of 
the experiment T = 0.1 K, there is no sign of a mag- 
netic transition in agreement with former ^SPJf^ and ac- 
susceptibility^^ measurements. The T-dependence of the 
total dc-magnetization M has been also measured be- 
tween 0.1 K and 3 K for fixed applied fields H in the 
range 0.05 - 0.7 T. Characteristic plots of M/ H versus 
T are presented in Fig. [21 At low temperature, satura- 
tion effects are evidenced by the decrease of M/H with 
increasing fields. More precisely, the field dependence of 
the magnetization measured at 0.2 K is plotted in the 
inset of fig. O At this temperature, the data are well 
described by the linear M = x(0-2 K)H relation for low 
fields < 0.1 T while saturation effects are clearly ob- 
served for H > 0.2 T. Therefore, the flattening at low 
temperature of M/H ~ x measured for H — 0.1 T in 
Figfflcan not be ascribed to a fleld effect. 

The NMR lineshift of chlorine in Herbertsmithite was 
measured recently^'"'. This local probe investigation is 



thought to give the intrinsic susceptibility Xi which is 
found to strongly deviate from the macroscopic SQUID 
data X for r < 50 K. Indeed, despite large error bars due 
to the broadening of the NMR line, Xi shows a broad 
maximum or at least a saturation below ~ 50 K. One 
can then put an upper limit to the intrinsic susceptibil- 
ity Xi < 1-5 X 10^'^ cm'^/mol Cu as represented by the 
shaded area in Fig. [TJ The maximum of Xi likely reflects 
a moderate enhancement of the short range AF correla- 
tions as in the well studied kagome bilayer casei^ii^i^S. 
Once these correlations have developed, it is doubtful 
that there will be a subsequent rise of the susceptibil- 
ity at lower temperature and we assume that the above 
upper limit for Xi ^iso stands down to K. The low 
temperature dc-susceptibility which is the subject of this 
report is therefore mainly dominated by a defect contri- 
bution. In the following we will make the simplest as- 
sumption that the intrinsic and defect contributions are 
uncorrelated and therefore X — Xd + Xi with Xi/Xd < 0-1 
for T < 2 K. As previously mentioned, one can antici- 
pate two types of defects in the structure which can both 
show a paramagnetic-like behavior and which both con- 
tribute to Xd in our analysis. First some Cu^"*" ions could 
lie on the interplane site. Their coupling to the kagome 
planes is likely very weak, maybe slightly ferromagnetic, 
as discussed in Ref^S'. Second, the dilution of the kagome 
magnetic lattice by Zn^+ ions is believed to locally stabi- 
lize dimers and to induce a weak staggered magnetization 
on further neighboring sites^^. This non trivial extended 
response of the system around a spin vacancy constitutes 
the second magnetic defect. In the closely related cop- 
per based anisotropic kagome structure of Volborthite^^, 
the controlled magnetic dilution by Zn/Cu substitution 
indeed yields a Curie-like tail that scales with the Zn 
contenbSS. 

We flrst consider the intermediate temperature range 
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1.5 K -10 K. As shown in the inset of Figll] 1/x — l/Xd 
does not extrapolate to when T — > as would be ex- 
pected for free spins following a Curie law. Instead Xd 
rather shows a Curie- Weiss behavior Xd = Cd/{T + 6d). 
A proper fit of the low T data requires an accurate knowl- 
edge of Xi{T)- In the absence of such data, we fit with 
a constant Xi in the two extreme cases; = which 
yields Od = 0.85 K and Q = 0.040 cmVmol Cu/K (fit 
range 1.5 K< T < 3 K) and Xi = 1-5 x 10"^ cm^/mol 
Cu which yields Od = 0.80 K and Cd = 0.0345 cm^/mol 
Cu/K (fit range 2 K< T < 10 K). A fit of the high tem- 
perature (T > 150 K) data gives a Curie- Weiss constant 
Ccw — 0.5 cm'^/mol Cu. If one assumes that the mag- 
netic defects behave as S=l/2 spins, their contribution 
corresponds to ^ 7% of weakly coupled S=l/2 spins out 
of the total Cu^^ contribution. This number is remark- 
ably similar to the estimated number of two level systems 
which contribute to the Schottky anomaly in heat ca- 
pacity measurementsiS^ and also of misplaced Cu^^ from 
neutron diffraction refinemen t This latter finding 
suggests that the main contribution to Xd comes from 
the interplane Cu^+ (S=l/2 defects) whereas the inte- 
grated staggered magnetization around a Zn^+ amounts 
to a rather small moment. 

More puzzling is the behavior below 1 K where a subse- 
quent enhancement of the susceptibility appears (Fig. [U 
inset). The above described Curie- Weiss regime accounts 
then only qualitatively for the flattening of x(T) (see 
dashed line in main panel and inset). At 0.1 K the rise of 
X with respect to the extrapolated Curie- Weiss behavior 
is about 1.3 x 10^^ cm'^/mol Cu, i.e. one order of mag- 
nitude larger than the upper limit of Xi- This enhance- 
ment is therefore also related to the defect contribution 
Xd- No Field Cooling-Zero Field Cooling opening could 
be detected below 1 K. Moreover x does not show any 
peak or divergence that would signal long range order- 
ing. Thus, the rise of x foi' T ~ 9d, probably reflects 
a strengthening of some ferromagnetic-like correlations 
between the magnetic defects rather than some kind of 
ordering. It is noticeable that a slight slowing down of 
the electronic spin fluctuation is detected at this same 
temperature T ~ in /iSR experiments. This nicely cor- 
roborates the correlation strengthening picture. Od is also 
close to the temperature of the maximum of the Schot- 
tky anomaly in zero field heat capacity data. Therefore, 
ksOd appears as a new energy scale for Herbertsmithite, 
most likely related to the magnetic defect system. 

The magnetic response of the defects strongly domi- 
nates the total susceptibility at low temperature and it 
is difficult to extract any information on the KAH con- 
tribution. However, one can expect different field de- 
pendences for the two contributions. Namely, the weakly 
coupled magnetic defects should be more easily saturated 
than the Cu^"'' spins belonging to the perfect kagome 
network with the stronger J ~ 190 K coupling. To fur- 
ther investigate the field dependence of the magnetiza- 
tion we measured M{H) curves in a vibrating sample 
magnetometer up to 14 T for constant temperatures in 
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FIG. 3: (Color online) Normalized magnetization {Msat = 
5583 emu) measured versus field in a VSM for 3 character- 
istic temperatures. The magnetization data at 0.2 K are 
also reported from the inset of Fig. [5] In the inset, for 
the same temperatures, M — XiH versus H/ (T + 9) with 
Xi = 1.25 X 10"^ cmVmol Cu and 9 = 1.3 K. 
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FIG. 4: (Color online) Full squares : normalized magnetiza- 
tion of Herbertsmithite measured at 1.7 K versus field. The 
solid line is a linear fit of the data for // > 10 T which likely 
reflects the non saturated intrinsic susceptibility Xi- Open 
squares : the defect contribution obtained by subtracting the 
above linear contribution from the full square magnetization 
curve. Dashed line : Brillouin function for 7.7% of free S=l/2 
spins. Inset: the shaded area on the SQUID x(^) plot repre- 
sents the possible values of Xi from this study. 



the range 1.7 K-25 K. Characteristic results are shown 
in Fig. [31 At 1.7 K, a strong saturation effect is ob- 
served above ~2 T and up to ^8 T where M{H) reaches 
a linear regime. At higher temperature, the saturation 
effect gradually disappears and at 10 K a nearly linear 
dependence is recovered. This behavior is compatible 
with a simple decomposition of the magnetization into a 
defect and an intrinsic contribution M = Md + Mi. One 
then assumes that for moderately high fields {H < 14 T) 
with respect to the coupling energy scale J, the linear 
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M,{H,T) = Xt{T)H relation holds and that the defect 
magnetization follows a Brillouin like saturation. In 
this scenario, at 1.7 K, the regime for iJ > 10 T where the 
magnetization is a linear function of i7, is explained by 
the complete saturation of the magnetic defects and the 
slope of M{H) is a direct measure of the intrinsic suscep- 
tibility Xi(1.7 K). In Fig.|3]the straight line corresponds 
to Xi(1.7 K) = 1.5 X 10~^ cm'^/mol Cu. It is the largest 
possible value for the intrinsic susceptibility so that the 
remaining defect magnetization extracted from our data 
(open squares) does not decrease at high fields. It is 
noticeable that the fully saturated defect magnetization 
amounts then to ~ 8% of the saturated magnetization 
of one Cu^"*" mole, in perfect agreement with the ~ 7% 
estimate given by the low T Curie- Weiss behavior of the 
defect susceptibility. However, one cannot exclude that 
the complex magnetic defects at play in Herbertsmithite 
are not completely saturated even at the lowest temper- 
ature and highest field of this study. Part of or the whole 
linear regime could then be ascribed to the defect con- 
tribution. The extracted Xi value is therefore only the 
upper limit of the kagome susceptibility at low tempera- 
ture. The possible values of Xi for T > 1.7 K from this 
analysis are represented by the shaded area in the inset 
of Fig. H 

As shown by the dashed line in Fig. 01 a simple S=l/2 
Brillouin function fails to capture the field dependence 
of the defect magnetization Md ~ M — XiH. Possi- 
ble reasons for this are that 1) the magnetic defects are 
complex objects involving the point defect itself, likely 
a misplaced Zn/Cu atom, and the local screening of the 
defect by the neighboring spins, so that one does not 
expect a simple S—1/2 effective spin value, 2) the mag- 
netic defects are slightly antiferromagnetically coupled 
which tends to reduce the field effect with respect to the 
H/T dependence of free spins. Note that unconstraining 
the spin value of the Brillouin function does not give 
either a good fit of the H and T dependence of Add- 
As shown in the inset of Fig. [31 the Md{T,H) data for 



1.7K<T<10K merge on a same curve if one uses the 
scaling variable H/ {T+6) which accounts phenomenolog- 
ically for the AF coupling. Good scaling is obtained for 
Xt = 0.00125 ± 0.00025 cm^/mol Cu and 6* = 1.1 ± 0.2 K 
in agreement with Ocw extracted from the low T Curie 
Weiss fit of the susceptibility. It is noticeable that be- 
low 1 K, deviations from this scaling appear gradually. 
Eventually, the 0.2 K curve can not be made to fall on 
the T > 1 K ones, even with different Xi and 9 values. 
It suggests that the effective defect moment does change 
below 1 K which corroborates the enhanced correlations 
scenario drawn from the analysis of the T dependence of 
the susceptibility. 

In summary, from a detailed study of the temperature 
and field dependence of the magnetization at low T , we 
can draw a coherent picture of the Herbertsmithite mag- 
netic behavior. The Curie-like tail in the susceptibility 
can be safely attributed to a defect contribution. The 
magnetic defects, probably of two kinds, behave in aver- 
age as weakly coupled spins {S ^ 1/2)- Signature of the 
coupling energy ~ 1 K are found ubiquitously in thermo- 
dynamics measurements as well as in the spin dynamics. 
The complex nature of the defects challenges both chem- 
istry to achieve a better control of Zn/Cu site occupation 
and theory to describe their magnetic behavior. Remark- 
ably, such a large quantity of defects does not seem to 
alter the underlying KAH physics. Besides, although the 
effect of the ~ 10 T external fields used in this study 
is not clearly known, our results are compatible with a 
finite kagome susceptibility at T ~ J/100 and thus ques- 
tion the ground state nature and the presence of a gap. 
Although we showed in this study that the low T up 
turn of the macroscopic susceptibility can be explained 
in a defect scenario without Dzyaloshinsky-Moriya per- 
turbation terms contrary to the initial proposal of RefJ^, 
they could nonetheless impact the low T intrinsic proper- 
ties and possibly increase the polarisability of the ground 
state of this unique realization of a S=l/2 kagome sys- 
tem. 
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